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Abstract 

We prove the well-posedness of the Cauchy problem on torus to an 
eletromagnetoelastic system. The physical model consists of three cou- 
pled partial differential equations, one of them is a hyperbolic equation 
describing the elastic medium and two other ones form a parabolic sys- 
tem, which comes from Maxwell's equations. Experimental measurements 
suggest that the elastic medium has a periodic structure, moreover with fi- 
nite number of discontinuities on the fundamental domain. Thus we have 
study in this paper the problem which we have defined as periodically 
Cauchy diffraction problem. 
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1 Introduction 



In this paper we provide a general framework leading to electromagnetoelastcity 
theory. In particular, this general theory encompass the two most important 
models given by a quasilinear and a semi-linear system of partial differential 
equations, described respectively by equations (see below) (|2.ip . (|2.3[) . (|2.4[) , 
(EH), dHU) and, (JOJ, ([231) . (EU), Moreover, considering the semi- 

linear case, and standing to plane waves, we have proved existence of weak 
solution, uniqueness and stability, therefore established the well-posedness of 
the periodically Cauchy diffraction problem, see Section l3?2l In fact, it seems 
that the notion of periodically Cauchy diffraction problem is one of the most 
realistic one as geophysical experiments suggest, that is, the medium has a 
natural periodicity with a finite number of discontinuities on the fundamental 
domain. Then, it is more physical correct to assume that the medium is periodic 
in space. 

The mathematical electromagnetoelasticity theory describes the interacting 
effects of an elastic solid medium and an electromagnetic field applied on it. 
Thus deformations experimented by the elastic solid are due to external electro- 
magnetic forces. More precisely, if an elastic electroconductive medium is imbed- 
ded in an electromagnetic field, then the elastic waves propagating through the 
medium will excite oscillations of the electromagnetic field and themselves will 
change under influence of the latter. Moreover, the waves which arise as a re- 
sult of such an interaction are called as electromagnetoelastic waves. We stress 
that, the first attempts to apply the theory of electromagnetoelasticity to the 
investigation of the wave propagation process in electroconductive media were 
made by Knopoff [12], Chadwick 4 , Dunkin and Eringen [7]. 

In fact, because of the importance of the applications, in particular to geo- 
physics applied in seismology, the petroleum reservoir research, the theory of 
electromagnetoelasticity have been developed fast recently. We address some of 
the important mathematical works on the propagation of electromagneticelastic 
waves, those are: Avdeev, Goryunov, Soboleva, and Priimenko pQ, Lorenzi and 
Priimenko |15) , Lorenzi and Romanov |16j , Priimenko and Vishnevskii [201 121] , 
Romanov [22l [23]. 



2 Non-linear electromagnetoelasticity theory 

The models of electromagnetism and elasticity could be given in a rational 
continuum physics way through the respective stored energy functions. First, 
let us see the electromagnetism theory, where we follow the Coleman and Dill 
work, see [5]. Let (t, x) g IxK 3 be the time-space domain and, we consider the 
magnetic intensity field H , the electric intensity field E, the magnetic induction 
B, the electric induction D, all of them taking values in M 3 . Moreover, we 
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consider the electromagnetic store energy function tp(D, B), thus we have 

d t D — cm\ x (dB*P) = —J, 

d t B + cwl x (d D iP) = 0, (2.1) 
div^ D = p e ,div x B = 0, 

where J is the induced electric current density field, also taking values in M 3 , 
and p e is the electric charge density. Further, we obtain from the stored energy 
function the following constitutive relations 

d D J(D, B) =: E h d B MD, B) =: H t (i = 1, 2, 3). 

We recall that, the first and second equations in (I2.1[) are respectively the Am- 
pere and Faraday's Law, and the last ones are constrains, which are compat- 
ible with the first two ones. In the linear theory of electromagnetism, called 
Maxwell's equations, we have 

D = eE and B = p e H, (2.2) 

where e, p e are respectively the dielectric and permeability tensors. There are 
many reasons to avoid the linear case, for instance the singularity at the origin 
on the electric field E. Therefore, some non-linear models have been proposed 
and the most famous one is due to M. Born and L. Infeld, see [2]. In this model 
the electromagnetic stored energy function is given by 

i>{D, B) = ^l + \D\ 2 + \B\z + \DxB\ 2 . (2.3) 



Now, we turn our attention to the elastic non-homogeneous medium. Let 
U(t,x) be the elastic displacement field taking values in M 3 . For elastic mate- 
rials, the constitutive relation of the stress tensor is given by 

T = T(F), 

where F is the gradient deformation of the solid, satisfying det F > 0. Moreover, 
considering isothermal deformations, which means that, the temperature of the 
medium is assumed constant and uniform throughout the entire solid, the stress 
tensor T could be written as 

T(F)+^/ d = T(F)=p^^F T , (2.4) 

where it is a constant, p is the mass function per unit volume of the solid, and 
Vl 7 is the elastic stored energy function. Furthermore, we have assumed that the 
elastic material is incompressible. Moreover, considering that the medium is 
locally isotropic and standing for linear elasticity, we could write 

T= A tr(S)J<j + 2/xS, (2.5) 
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where A, p are scalar functions called the Lame elastic moduli and S is the 
infinitesimal strain tensor, given by 

S := i(VC/ + VJ7 T ). 

We recall that (|2.5p is also known as Hooke's law. Therefore, from the Cauchy's 
first law and the Hooke's law, we have the following equation describing the 
evolution of the linear elastic electrically-magnetic conduction solid 

pdf t U = div x 7+F b , (2.6) 

where F b := J x B + p e E + F s is the body field force per unit mass, more 
precisely J x B is the electric- magnetic part that retards the motion of the solid, 
p e E stands the part of the body force due to the existence of the charge density 
p e , and finally F s are others body forces, for instance gravitational effects on 
the solid. 

Therefore, we have seen that the electromagnetic field influences the elastic 
field by entering the elastic stress equations of motion as a body force called 
Lorentz's ponderomotive force. Now, we are going to establish that the elastic 
field turn influences in the electromagnetic field by modifying the Ohm's law. 
Indeed, in a moving conductor medium the current is determined by Ohm's law, 
see [12], that is 

J = a (E + d t U x B) + p e d t U + F e , (2.7) 

where a is the electrical conductivity tensor field and F e is an external elec- 
tromagnetic force. Then, we see that the current distribution is modified by 
the elastic deformations. Thus the interaction between the elastic field and 
electromagnetic field is expressed through equations (|2.6|) and (|2.7j) . 

Now, we are in position to establish the basis of the electromagnetoelasticity 
theory. In fact, equations (|2.ip. (12^21 . (|2.5p. ((231) and (j2~7) form a semi- linear 
system of partial differential equations and, (|2~Tj) . (|2l3l) . (|2~4l) . (|2J)) and (|2~7) a 
quasi-linear one, which is much more complicated, since shocks are allowed to 
exist, see Dafermos [6], and for instance, Neves & Serre [18] for the nonlinear 
Maxwell's equations setting. In this paper we consider the semi-linear case 
and in other to solve the problem both mechanical and electromagnetic data 
should be given. We stress that even in the semi-linear setting, the nonlinearity 
involved is non-trivial and difficult to hand up. 

Finally, we observe that the theories of electromagnetoelasticity have been 
extended in various ways. For more acquaintance with the modern state of the 
theory, the reader is addressed to, e.g., j3][8]. 

An outline of this paper follows. In the rest of this section we fix the notation 
and give some mathematical definitions, which will be used during the paper. In 
Section [3] we establish our problem, where we give the exact notion of Periodi- 
cally Cauchy diffraction problem. Finally, Section 2] presents the well-posedness 
result for the problem proposed at the preceding section. 
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2.1 Functional notation and background 

At this point we fix some functional spaces, which will be used in the paper. For 
convenience, we follow the notation of Ladyzhenskaia, Solonnikov and Uralceva, 
see [13] ■ Although, we will be in the periodic setting, which is not the case in 
[T3] . Therefore, we address also [TT] and [25] . 

We will concern periodic functions in R (i.e., in the spacial domain). For 
simplicity, we take the fundamental period to be one. So, we define T := R\Z to 
be the 1-dimensional torus, it means that, a scalar periodic function / : T — > R, 
satisfies 



Moreover, we recall that a periodic function is completely determined by its 
values in the fundamental domain, here we take ft := [0, 1). Therefore, periodic 
functions on R will be considered as functions on T or functions on f2. In fact, 
the point of view depends only on the context. 

By dz, we denote the Lebesgue measure on R and for measure-theoretic 
purposes, we do not distinguish between the measure induced on T and the 
Lebesgue measure, that is, if / is any function on T, we define 



when the right hand side makes sense; indeed a function / on T is measur- 
able when the corresponding function on ft is Lebesgue measurable. The same 
definition could be obtained with the concept of periodic distributions and a 
periodic partition of unity, see Vo-Khac Khoan [IT], Tome 2. 

The Banach space L q (T) consists of all measurable functions on T that are 
gth-power (q > 1) summable on T provided with the norm 



For some T > 0, we set n T := (0, T) x T and analogously Q T := (0, T) x fl 
The Banach space L,j j7 (n<r), (q, 7 > 1), consists of all measurable functions on 
with a finite norm 



Moreover, when q = 7 the Banach space L q q (I\T) will be denoted by L g (IIx)i 
and the norm ||i>||g,g,n T " by |M|q,n r - 

Weak (generalized) derivatives should be understood in the customary way, 
we address the book of Evans, [9]. For / an integral and q > 1, we denote by: 



f(x + k) = f(x) for each i£l and k£Z. 
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W l q {T), the Banach space consisting of all functions of L q (T) having weak 
derivatives of all forms up to order I inclusively, that are gth-power summable 
on T. The norm in W'(T) is defined by the equality 



MI?i = Ell^ll,.T. 



s=0 

I 



W q (T), the closure in H^(T) of all functions that are infinitely differentiable 
and finite in T. 

W? ' 1 (JIt), the Banach space consisting of the L g (nT)-elements having weak 
derivatives of the form D r t D s z with any r, s satisfying the inequality 2r + s < 21. 
The norm in it is defined by the equality 

N&=E £ W d td>\kn T , 

j—0 2r+s—j 

where the summation ^ 2 r+s=j IS taken over all nonnegative integers r and s 
satisfying the condition 2r + s = j. 

W2 ' (n T ), (k = 0, 1), the Hilbert space endowed with the scalar product 



(u,v) w i,k {nT ^ = // (uv + u z v z + ku t v t )dzdt. 

J J TIt 

V2(Ht), the Banach space consisting of all W / 2 1,0 (riT)-elements having a finite 
norm 

\v\u T = vrai max ||u||2,T + |K||2.n T , 
te[o,T] 

where here and below 

.1/2 



v 2 ||2,n T = ( // (v z ) 2 dzdt) 



7 2 1,0 (n T ), the Banach space obtained by completing the set W 2 ' (ILt) in the 
norm of ^(IIt). 

2 (IIt)) the subset of those elements v £ V^ 1 '°(IIt), such that 

T-r r 

1 (v(z,t + t) - v(z,t)) 2 dzdt -> asr ->• . 



Remark 2.1. C a ' a ^ 2 (QT) is the set of all continuous functions in Q T satisfying 
Holder conditions in the spacial variable z with exponent a and in the time 
variable t with exponent a/2. Following Stein & Weiss f25jj , we remark that 
C(T) does not correspond to the class of continuous functions on O, but only 
to those functions which remain continuous on R when extended by periodicity. 
Therefore, C a ' a / 2 (JlT) should be understood in this sense. Although, the spaces 
Lq(Q), L q ^(Qt), etc. are defined in a similar way as above, where we replace 
respectively T by fl and Ht by Qt- 
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3 Statement of the problem 



3.1 Plane waves 

We are going to focus on plane waves, which depend only on the one scalar space 
variable and time. Let us consider X3, the such spacial coordinate. Therefore, 
all of the fields involved in equations (jlTTj) . (|2T2"j) . (|2"3)) . (|2l))) and (12^7) depend 
on (t, 0,0, x 3 ) = (t,x 3 ) variables. Moreover, the vector fields F e and F s have 
the following representations 

F* = (Ff,Fi,0), F s = (0,0, Fi), 

where Ff, F 2 and F£ are scalar functions and, from now on, z stands for the 
variable X3. 

In the case of diffusion electromagnetic processes the time derivative of D in 
the Ampere's Law is very small in comparison with the conduction current J 
and, and in such way could be dropped. Further, we assume that p = const, p e = 
const, a is a scalar function and the charge density p e = 0, that is, the media 
satisfies the quasi-neutrality condition. Then, considering these assumptions, we 
have from (pTTj) . (j2~2j) . (j2~5l) . (|2~6)) and (|2~7)) the following system of semi-linear 
partial differential equations 

d t Hi = d z (—d z Hi -HxdtU- — F?), (3.8) 
d t H 2 = d z (—d z H 2 -H 2 d t U- — F£), (3.9) 

(T fi e Ope 

dlU = d z (v 2 d z U - ^{Hf + Hi)) + -Fi, (3.10) 
2p p 

where v := \/(A + 2p)/p is the longitudinal elastic wave velocity and U := U3. 
For simplicity of notation and similarity purposes, we will form a dimensionless 
system of equations. Let L, Vq and Ho be the characteristic values of length, 
velocity and magnetic field, respectively. We set, 

1 He Hi v 



PeLVoa' P: 2pV 2 ' " ' V ' 

respectively the first and third one, the dimensionless magnetic viscosity and 
the dimensionless velocity of the elastic waves propagation. Then, after simple 
transformations, we obtain from (|3.8[) - fl3.10[l 

h t = (rh z - hu t - rj) z , (3.11) 
u tt = {v 2 u z ~ph 2 ) z + f, (3.12) 

where h := (hi,h 2 )^ j := (ji,j 2 ), u and / are respectively the dimensionless 
analogues of 

H t ,F^ (i = 1,2), U, -Fi, 

P 
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and h 2 = h - h = h\ + 



3.2 Periodically Cauchy diffraction problem 

In this section, we state our initial-value problem on torus, that is, we con- 
sider equations (|3.11|) - (|3.12j) in IIt- We assume that p is a positive number, 
and r(z), v{z) are positive bounded 1-periodic functions, piecewise smooth and 
discontinuous at each point z = Zf. + K, (k = 1, 2, . . . , to; to € N), (k € Z), with 

< Zi < Z2 < • • • < 2 m < 1. 

The system (|3.11j) - (|3.12j) of partial differential equations is supplemented by 
the initial-data 

h = h on {0} x T, (3.13) 
u = uo, u t = ui on {0} x T, (3-14) 

where hy, uq and u\ are given 1-periodic smooth functions. 

Now, since the coefficients in (|3.11l) - (|3.12[) have to discontinuities in J7, we are 
going to say that (|3.11[l - (|3.14[) form a periodically Cauchy diffraction problem, 
which is defined in analogy with one at Ladyzhenskaia, Solonnikov and Uralceva, 
see p21 pp. 224-232]. Thus for each k — 0,1,..., to, we set fi fc := [zf.,Zf.+i), 
with zq = and z m+ \ = 1. Therefore, we have 

m m 

n = (J n k and Q T = U Qt, 

where Q\ = (0, T) x It follows that, in each domain Q\ there is given a 
parabolic-hyperbolic system with smooth coefficients and free terms. One of the 
main purposes is to find in ITr a (weak) solution (h, u) of this system, satisfying: 

• in Qj., k = 1,2, ... ,m, the corresponding equations (|3.11p — (|3.12|) ; 

• on the lower base of 11^ the initial condition (|3 . 1 3[) (|3. 14[) ; 

• at the jump points z = Z). + K, the following compatibility conditions, 

[h] = 0, [u] = 0, (3.15) 
[r(h z -j)}=0, [v 2 u z ]=0, (3.16) 

where [v] denotes the jump of the function v as it passes through a discontinuous 
point. 
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Remark 3.1. Here, we observe our strategy to show existence of solution to 
periodically Cauchy diffraction problem (|3 . 1 1 1) (|3 . 1 6[) . First, we are going to 
consider the associated initial-boundary value problem given by the equations 
(|3.1ip - (|3.12l) posed in Qt, the initial-data (| 3 . 1 3 1) (| 3 . 1 4 [) on {0}xf2 and the jump 
conditions ()3.15|) — (|3.16|) . Moreover supplemented with the following boundary 
conditions 

h(t,0) = h(t,l), M*»o) = MM), (3-17) 
u(t,0) = «(t,l), u,(t,0) = u a (i,l). (3.18) 

Once we have constructed the existence of solution to (| 3 . 1 1 1) — (| 3 . 1 8 1) in T)'(Qt), 
we extent it to Tlx by periodicity and show that it is a solution to the Cauchy 
problem (j3.11[) - fl3.16[) in Ilr, as the definition given below. 

The following definition tells us in which sense a pair (h, u) is a weak solution 
to (I33TD-(I3~TC1). 

Definition 3.2. A pair of 1-periodic functions in the z-spacial variable 

hG7 2 (n T ), ue^'prf 



is called a weak solution of the initial-value problem (|3.11|l - (|3.16|l if satisfies the 
identities 



hi rjt dzdt + rhi z r] z dzdt —II hi Utf]z dzdt 

1. 1 -/ J Jn T J Jut 

rjii] z dzdt+ / hi i](0)dz, (1 = 1,2), 



Ut Ct dzdt + / / ^ 2 u z Cz dzcft + / / P (h 2 )z C ckoW 
f(dzdt+ uiC(0)dz, u(0)=u , 



(3.19) 



(3.20) 



/or a/Z 77, C € VF 2 ' (Ilr), which are equal to zero for t = T. 

Moreover, it is also possible to define the weak solution of (|3.1ip - (|3.16p in a 
way somewhat differently. 

Definition 3.3. A pair of 1-periodic functions in the z-spacial variable 

f heVt(n. T ), uew 2 M (n T )) 
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is called a weak solution of the initial-value problem (|3.1ip - (|3.16p if satisfies, 
for almost all ti £ [0,T], the identities 

hi i] t dzdt + I / rhi z r\ z dzdt hiUtT] z dzdt 

(3 21) 

rjirj z dzdt+ / h w n(0)dz- / h t (ti) ry(ii) dz, (1 = 1,2), 
o Jt Jt Jt 



u t Ct dzdt + / / v 2 u z C, z dzdt + / p (h 2 ) z Q dz dt 

10 JT JO JT JO JT 



o Jt 



(3.22) 



f(dzdt+ / ui((Q)dz- / u t (h) ((h) dz, u(0) = u , 



for any -q, ( G W^'^IIt). 

Both definitions are equivalent. The fulfillment of transmission conditions 
(13.161) is understand in the sense of the identities considered in Definition 13.21 



4 Well-posedness 

The main focus of this section is to establish in the sense of Hadamard the well- 
posedness of the periodically Cauchy diffraction problem (|3. 1 ip — (13. 16p . that 
is, we prove existence of a solution, show the uniqueness and the continuous 
dependence in all the parameters and data. 

4.1 Existence of weak solution 

In this section we show a week solution of the periodically Cauchy diffraction 
problem (|3~TT|) - (|3T6| . 

Theorem 4.1. (Existence theorem). Let h e C a (T), (a e (0,1)) and 
(uq,ui) £ W2(T) x ^2(T) be given initial-data to the Cauchy problem (|3.11[) 
(|3.16p . Moreover, assume that the constant p, the functions r,v and the free 
members j, /, satisfy the properties: 

(a) r,v,j,f are 1-periodic functions, piecewise smooth, bounded and dis- 
continuous at each point z — Zk + k, (k = 1, 2, . . . , m; m £ N), (k £ Z), 
with 

< Zx < z 2 < ■ ■ ■ < z rn < 1; 

(b) p is a positive number and there exist < To, Vq, T\, V\ < oo, such that 

ro < r < r\ and Vq < v < V\. 
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Then, the Cauchy problem (|3.1ip - (|3. 161) has a weak solution 
h€V 2 (Il T ), ueW^QIt) 



Proof. 1. First, following the strategy given at Remark 13. H we consider the 
system given by (|3.11[) - (|3.12j) in V'(Qt), the initial-boundary value data given 
respectively by (|3.13P — (|3.14p and (I3.17p - Q3.18p and finally the jump conditions 
(|3.15p - (|3.16p . To solve this initial-boundary value diffraction problem we make 
use of Faedo-Galerkin's method. 

2. Now, we consider in W^i^l) a fundamental system of functions ipk{z), 
(k = 1, 2, . . .), satisfying for each k > 1 

^(0) = ^ fc (l), fe(0) = i>kz(l)- (4.23) 

Moreover, we assume that {tpk}kLi is an orthonormal basis of La(n) and an 
orthogonal basis of W^ity- Fixed a positive integer N, we will look for an 
approximating solution in the form 

N N 

h N (t):="£a%(t)i; k , u N (t) := £ &f (i) (4.24) 
fc=i fc=i 

where for < i < T. The coefficients aj^ = (a^,, a^,), 

«»(*) = M bg(t) = (u N (t), tpk) ,1 = 1,2 ,k = l,. ..,N , 

are determined from the equations 

(fcf , = ( - rfc£ + /if uf + rj u i> kz ) , 

{h?(Q),ip k ) = h m ,l = l,2,k = l,...,N, 



(4.25) 



(u N ^ k )tt = (~v 2 u? + P h 2 ,4, kz ) + (/, Vfe), 
(u N (0),^ k ) = M fe, (^(0), il> k ) t =u lk ,k = l,...,N. 



(4.26) 



Moreover, the values ho k = (hio k ,h2ok),uo k ,ui k are the Fourier coefficients in 
2^2 (fi) °f initial data with respect to the system of the functions ip k . 

3. The equations (I4.25[) and (|4.26p form a system of nonlinear ordinary 
differential equations, and from standard existence theory, there exists a unique 
local solution on a maximal interval [0, r), for some r > 0. In fact, we will prove 
that |dfc |, I; & — 1> • • • > 2V, are bounded functions for £ 6 [0, T]. Therefore, the 
system (|4.25l) - (|4.26p has a unique solution on [0,T] for any T > 0. Indeed, we 
multiply the differential equations in (|4.25[) and (|4.26[) respectively by p af k and 
6^. Now, we sum the obtained equalities over all k from 1 to N and integrate 
the result with respect to t from to t\. Therefore, we obtain 



^{plKmln + \\u?(t)\\ln + \\™?(t)\\ln) "_ n + p\\^ \\l Qti 



— AT II 2 



t=0 



2 

"I, 



= P I I rj-h N z dzdt + /wf dzdt, (4.27) 
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where Qt 1 = (0,tx) x il. One observes that 



A' 



l|/A0)|li,n = £«I(0)<IIMi,n> 



fc=i 



2 



ll^(0)|||n = £&L(0)<N| 
fc=i 

lk«f(0)||2,n <W)"?ll«o.|la,n, 



ri 



where the positive constant /io does not depend on N. Consequently, it follows 
that 

fl|fc"(*i)ll!,n + ~lk"(*i)lll,n + 5ll^(*i)Hl, n + ~\\Vth?\\l Qti 

,N 



where the positive constant 



<mi+ / / \M\dzdt, 

0*1 



Pill. 112 , 1 IU. ||2 , Win , |,2 P / / . -2 



/*i = |ll^llln+2ll u illln + ^ll«o»||2,n + | / / rj'dwft 
does not depend on N. In particular, we have 



< Mi + I // /<dsdt| (4.28) 

Otj 



2 



and integrating this inequality with respect to t\ from to T, we obtain 



T <»iT+ [ \ff fu? dzdt\dh. (4.29) 

JO JJQ tl 



\\\<wi 



Applying the generalized Young's inequality, we have 

fu?dzdt\<^ [[ (u?) 2 dzdt+^- [[ fdzdt (4.30) 

and setting e = 1/2T, we obtain 

/ | f [ /uf dzdt\ dh<~ [[ (uf ) 2 dzdt + T 2 [[ f dzdt . (4.31) 



Therefore, inequalities P~2"9")l and (|43T|) yield 
1 



4ll«t llio* <HiT + T* If fdzdt. 
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From the latter inequality and (|4.30p with e = 1/2T, it follows that 



<Hi+2T 



f dzdt . 



Consequently, from (|4.27[) we have for any tj € (0,T] 



ll^(ii)ll^+2lK(*0III,n + 9lk<(*i)III,n+Pll^llI,Q tl 



<2^ + 2T I f z dzdt =: fi 2 , (4.32) 



where the positive constant \ii does not depend on N. It follows from (|4.32|) 
that , 6^, 6^ are uniformly bounded functions for each t £ [0,ii], t\ < T. 

Remark 4.2. One observes that, for each t £ [0,ti], t\<T 



h N (t,z) - h" (t))dz = 0, 



where h N is the average of h N on ft, that is 



h N (t) := / h N {t, z) dz < \\h N (t)\U,a < ( —) 
Jn v P ' 



1/2 



4. Now, let us show that for an arbitrary fixed k < N, the a^,b^,b^ t are 
also equicontinuous on [0,T]. For any S > 0, we set Qt,t+s = (t,t + 5) x fi. 
Then, we have from (I4.25[) 



a%(t + 6)-a%(t)\ = 



N i 



(-rh™ + u?h N + rj) ^ kz dzdt 



< 



+ 



If \rh^ kz \dzdt+ ff |< h N \p kz \dzdt 



' Qt,t+s 

\rj ip kz \dzdt =: h + 1 2 + h, 



t,t+<5 



with obvious notation. For convenience, we recall generalized Holder's inequality 

3 

\viv 2 v 3 \ dzdt < |IMsj,7j,Qt> ( 4 - 33 ) 



3=1 



where 



111 111 

q u 7l £ l,oo), i = 1,2,3, — + — + — = 1, — + — + — = 1. 

qi <?2 <73 71 72 73 
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Then, applying the above inequality, we could write 

h < n\\h* h,Q ht+s ||^fca||2,Q t , t+4 
\pr J 

h < ||2,4,Q t , t+s \\h N \\oo,4,Q t:t+s \\lpkz\\2,Qt, t+s 

<M2 S 1 / 4 II^IUaq^.^ <S 1/2 ||lMa,n 
< [ii S 3/4 \\h N \\ 

^3 < ri||j|| 2 ,Q t)t+a \\^kz\\2,Q t , t+s 

<n [i 2 tf 1/2 ||Vk*lkn. 

In order to estimate \h ||oo,4,Q M+a , first we recall Remark [4T2l and define 

v := h — h . 

Then, we apply Theorem 2.2. in [T31 pp. 62-63], with Remark 2.1 in [T31 pp.63] 
for each component of the vector- function v . Therefore, we have 

a i_a 
||v|| g ,n < 0IKIl2,n ' INk^ 

where 111 
0€[2,+oo], 7 e[4,+oo], - + — = -, 

7 2g 4 

the positive constant ft depends on 7 and f2, and v stands for the components of 
v . By integration with respect to the time variable from t to t + 5, we obtain 

||«||g,7,g«.«+* < 0IMIa,Q t>t+i Vrai re ^ <5] H U ll2,f2 T - 

Moreover, using Young's inequality, we could rewrite the latter inequality as 

2 2 

IMI<r,7,Q t ,*+i ^ 0~IKIkQ t , t+6 + /3(1 - -) vrai max ||«|| 2 ,n . (4.34) 

7 7 re[t,t+oJ 

Consequently, from equation (|4.34[) . it follows that 

,JV,,l/2 . . _ 



/ 2 < M 2^ /4 /3(||^||^, t+5 +vrai^^ ||^|| a , n 
+ /i 2 ^ 3/4 /? H^Qt.t+i + V1 ' ai max II^MII^J HV^Ikn 



<2M2^ 4 /3((^) 1/4 + (^) 1/4 )||fe||2,o. 
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Therefore, given e > 0, there exists a S > such that, if0<i<r<t + <5, 
then 

|af(r)-af(«)|<e||^|| 2i n, 

where e does not depend on N and tends to zero as S — > + , i.e. the uniformly 
equicontinuity in t of the . The equicontinuity of functions follows from 
the boundedness of their derivatives. Moreover, it is proved similarly as done for 
the functions , N > k, that b^ t are equicontinuous functions on t for N > k. 

5. By usual Cantor's diagonal process, we can select subsequences 
{a^}£ =1 c{af}£ =1 and {6f"}£ =1 C {&f }£U 



be converging uniformly on [0,T] to some continuous functions a,k(t),bk(t). 
Hence for each t £ [0, T], we define 

oo oo 

fc(t) := <**(*) V>fc, u(t):=5^6*(*)^ k . (4.35) 
fe=i fe=i 

The sequence {/i^" 1 } converges to the function h weakly in L2^l) and uniformly 
with respect to t in [0,T]. Indeed, for any ijj £ L2(&), taking ip = ipe, with 
e = (1,1), it follows that 

s OO 

= J2^^k)(h N ™-h,ip k e) + ( Y, (h N ™-h,(iP,i> k )i> k e), (4.36) 



k=l k=s+l 



with 



\(h Nm -h, J2 (lMkW*e)| <Cx( ^ faifitf) 1 '* = 01.11(8), 

k=s+l k=s+l 

where the positive constant Ci does not depend on N m and s. Now, we choose 
s large enough so that, CiR(s) becomes less than a preassigned 2e > 0. On the 
other hand, for fixed s and large enough N m , the first sum in (|4.36|) will be less 
than e for all t in [0, T]. Therefore, the term 

\(h N ™-h,^)\ 

can be made less than e for all t in [0,T]. Consequently, the sequence {h N "'} 
converge to h weakly in L^^l), uniformly with respect to t £ [0,T]. Moreover, 
by construction and uniqueness of the limit, we obtain for each t £ [0, T] 

h(t,0)=h(t,l) h z (t,0) =h z (t,l). (4.37) 

The sequence {u Nm } is bounded in Loo(0, T; W / 2 1 (17)) and the sequence {u^ m } 
is bounded in Loo(0, T; L 2 (^))- Hence 



converges to u *-weakly in L oo (0,T; W 2 1 (^)), 
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u t m converges to u t *-weakly in Loo(0, T; L 2 (f2)). 

Moreover, the sequence {/"■} is bounded in W / 2 1,1 (Qt)- Then, from the Sobolev's 
Imbedding Theorem the sequence {u Nm } converges strongly to u in L2(Qt) and, 
passing to an appropriate subsequence if necessary, we obtain a.e. convergence 
in Qt- Finally, we observe that for a.e. t £ [0, T] 

u(t, 0) = u(t, 1) u z (t,0) = u z (t,l). (4.38) 

6. At this point, we proceed to extend by periodicity the considered func- 
tions, where conditions (j4.37[) and (|4.38p will be used. For convenience, we 
remain the same notation. Therefore, we have 

h Nm (t) converges to the function hit) weakly in Z 2 (T), (4.39) 

uniformly with respect to t in [0, T] and 

/- converges to u *-wcakly in £«,((), T; W^{T)), (4.40) 

u^ m converges to u t *-weakly in L^O, T; L 2 (T)). (4-41) 

Now, let us show that the sequence {uf } is bounded in L 2 (0, T; iT _1 (T)), where 
iJ-^T) is the dual space of H 1 ^) = W 2 1 (T), see Iorio-Iorio [TO]. For this 
purpose, we fix any function * £ W%(T), such that ||*||wri(T) < 1, * = *i +*2, 
where 

*i € spanjVfcljfcf , and (* 2 , fa) = 0, k = 1, . . . , N. 

We denote by (a, b) the pairing between a £ _ff _1 (T) and b £ W^T). From 
(|4.26j> . we have 

(u£, *) = «, * x ) = -(^f , * u ) + p((h N )\ 9 ls ) + (/, tt x ) . 

It is easy to check that H^ilj wi(T) < 1- In view of this and estimates obtained 
before, we obtain 

|<u£,tt>| <C 2 /i a , 

where the positive constant C2 docs not depend on N. It proves that {u$} is 
bounded in L 2 (0, T; i/ _1 (T)). Then, it follows that 

Uu u tt weakly in L 2 (0, T; fl" _1 (T)) . (4.42) 

Analogously, let us show that the sequence {/if} is bounded in Z, 2 (0, T; ff-^T)). 
First from (|4.32p , one can extract from the sequence { h m } a subsequence con- 
verging to h weakly in L2(IIt) together with {h^ m }. Now consider a function 
$ £ Wi(T), such that H^llwjpD < 1, * = $x + $2, where $1 £ span{#}^f , 
and (&2,i>k) = Q,k = 1,. . . ,N. Since we have 

AT 

&tt = E a ^ fc € span{^}fe:f , 1 = 1,2, 
fc=i 
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hence the following equalities are valid 




rji,®i z ), 1=1,2. 



Using ||$i|| w i(n) < 1 and l[02"]) - (|£3"3p . we obtain 

K^,$)|<C3M2, ^ = 1,2, 

where the positive constant C3 does not depend on N. The latter provides that 
{/if} is bounded in L 2 (0, T; ii" _1 (T)). Also from (IQ21) . it follows that {h N } is 
bounded in £ 2 (0,T; W^T)). The well-known Aubin-Lions Theorem, ;14, p.58], 
implies that, a subsequence of {h^} converges strongly in L 2 (IIt). Therefore, 
passing a subsequence if necessary, we have a.e. convergence in IIt- Without 
loss of generality, we can assume that the sequences {h m }, {u Nm } converge 
to the limit functions /i, it in the sense mentioned above. Consequently, the 
sequence {h N } converges a.e. in IIt. 

7. Finally, it remains to show that the pair (h, it) of limit functions satisfy 
equations (|3.19l) and (|3 . 201) . First, we will show that the function h satisfies 
equality (|3.19p . For this purpose we multiply each equation of (|4.25[) by a 
smooth function afc(t) that is equal to zero for t = T, then sum over all k from 
1 to N' < N, and integrate the result with respect to t from to T. Therefore, 
we obtain 



+ K,?f (0)),/ = l,2, (4.43) 
where T N '(t,z) = EkLi a k(t)M^) belongs to ^(0, T; ^ L^Ut). 



We claim that we can pass to the limit in equation (|4.43[) along the subsequences 
selected above, assuming T N fixed, and thereby arrive at (|4.43p with h^ m , u Nm 
being replaced respectively by hi,u, i.e. we must prove that 



But the first term in the right hand side of the latter equality tends to zero as 





as m 00 




Indeed, by a simple algebraic manipulation, we have 
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m — > oo according to the estimates 
Jo 

< nf iu TO (n T ) [ T \\u^(t)\\L 2m \\(h^ - h.mh^dt 

Jo 

< 2 M2 T||Tf \\ Loo {n T )\\(h? m - MIUoo(o,T;£ 2 (i)) 0, 

as rn — > oo. Moreover, one observes that the last term in (|4.43[) tends to 0, 
since we have proved that u^ m converges to u t *-weakly in Lqo(0, T; L 2 (T)) and 
h t rf e L 2 (0,T;L 2 (T)) ^ Li(0,T;L 2 (T)) due to the estimates 

\\l 2 (0,T;L 2 (T)) < \\f l l\\L 2 (.0,T;L oo (T))\\'Y^ II L oo (0,T;L 2 (T)) i £=1,2. 

But it is not difficult to show that T N are dense in the space of all required 
functions in the first definition of a weak solution. In view of this hi, I = 1,2, 
satisfy equalities (|3.19p and belong to be the space V^Hr)- Moreover, since 

Nflli,n<2/i2, 

passing to the limit, we obtain 

Klll.il <2M2- 

Therefore, we have 

u? e£ 1)TA _(II T ), K€(0,-), 

and we can deduce from Theorem 8.1 [131 p. 192], interpreting Q = [0, T] x [— 1, 2], 
Q' = [0,T] x [0, 1] and by the periodicity procedure, that 

vraimax < Ca, 
n T 

for some positive constant C4. Similarly, from Theorem 10.1 of [131 p. 204], it 
follows that 

where < u\ < a. Now, let us show that u satisfies equality (13.201) . For this 
purpose, first we multiply each equation of (|4.26[) by the ak(t), then sum over 
all k from 1 to N' < N, and integrate the result with respect to t from to T. 
Hence after an integration by parts, we have 

f{u» , Tf ) dt = f ((Su» , Tf ) - p((h N ) 2 , Tf ) + (/, T N ')) dt 
Jo Jo 

-(uf,T w '(-,0)), u N (z,0) = u»(z), (4.44) 

where T N were defined in (|4.43[) . Then, passing to the limit in (|4.44l) with 
respect to the subsequence {N m } selected above, assuming that T N is fixed, 
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we arrive at (|4.44[) with h N ,u N being replaced by h,u. Since {h N } converges 
a.e. in Ht to h, we can pass to limit in the non-linear term {(h N ) 2 }, which 
converges to (h) 2 weakly in (|4.44l) . Moreover, since maxn T \h\ < C, we have 
that 



(h 2 ,rf)dt} 



is bounded for any T N £ W^iJ^r)- Therefore, as the functions T N are dense 
in the space considered in the definition of weak solution, we conclude that 
the function u satisfies equality p. 201) and is a weak solution from W 2 ' (ITr)- 
Consequently, Theorem 14. II is proved. □ 

Membership of such a solution h in ^'^(Ht) follows from Lemma 4.1 of 
p~3] p. 158] and Theorem 14.11 Then, we have the following 

Corollary 4.3. Any weak solution h of problem (|3. 1 1[) — (13. 16|) from V2 (ILr) 

belongs to ^'^(IIt)- 

Moreover, we have also that 

Corollary 4.4. For any function <f> £ L 2 (0, T; W^CT)) is valid the following 
equality 

(u tt ,(f>)dt= / ((f -p(h 2 ) z ,4>)- {v 2 u z ,(f> z ))dt. 
Jo 

The latter result allows us to conclude that the following equation holds for 
any function £(z) € W^T) an d f° r almost all t € [0,T], that is 

{u tt ^) = {f-p{h\,0-{v 2 u z ^ z ). 

Furthermore, we notice that 

u e C([0,T];L 2 (T)), u t e C([0, T}; H^T)) . 



4.2 Uniqueness of weak solution 

In this section we prove uniqueness of weak solution of the initial- value problem 
(|3.11|) - (|3.16p . First, we need the following 

Lemma 4.5. Suppose that 

heV 2 (n T ), ueW^iUr) 

is a weak solution of problem (|3.11[) - (|3.16j) . Then, the following inequality is 
valid for almost all t\ G [0, T] 

i f (ph 2 (ti)+u l {t l ) + v' 1 ul{t l ))dz+)- f f prh 2 z dzdt 
* Jt * Jo Jt 

< I ' (ph 2 + u 2 1 + v 2 u 2 )z )dz + [ [ prj 2 dzdt + 2t 1 [ [ f 2 dzdt. (4.45) 
jt Jo Jt Jo Jt 
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Proof. 1. Let e G (0,T) be fixed and consider the following functions 
Ve(t, z ) = ~ v(t, z) dr, Ce(t, z) = - / C(r, z) dr, 



where 77, C € W 2 ' ((-e,T) x T), with 

r)(r, z), C(t, z) = for t G [-e, 0] U [T - e, T}. 
Then, we take respectively fje, Cc as V an d C m (|3.19[) (|3.20l) . Since 

(??e)t = (j)t)g, 

we could rewrite the first term in p,19p in the following manner 

T-e p pT-e r 

/ hi f/tt dzdt = — I huf/t dzdt 

1) Jt J ■ Jt 

T- 

hut f] dzdt, I = 1, 2, 







where 



1 



6 



t+e 



h e {t):=- h{T)dT, 1 = 1,2. 



Analogously, we have 



T-e [■ /-T-e r 

j u t C u dzdt = / u ett C dzdt. 

Jo Jt 



where 



1 



r 



t+r- 



u e (t) := - u(t) dr. 



2. Now, we proceed similarly in the remainder terms of p.l9|) - ()3.20p . to 
transfer the averaging operation (-) g from fjg, (,1 to respectively coefficients. 
Then, with the above notation and taking into account the permutability of the 
averaging operation with differentiation with respect to z, we obtain 



T- 



T-e 



(hut V + (rh lz -h t u t - rji) e jfej dzdt = 0, I = 1, 2, 

(4.46) 

(v 2 u ez ( z + (u tt + P(h 2 ) z - f) e C) dzdt = 0. 

Let ti G [0,T — e] be fixed. For convenience, we choose £ equals zero for each 
t > h and C = C on [0,ti] x T, for some C G W^CC ) *0 x T )- Analogously, we 
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take fj = for t > t%, and by standard density argument, f\ = r\ on [0,ii] x T, 
for some r\ € 

^((Mi) x T )- Therefore, we obtain from (j446|) 





7' 



JT 



{hut V + (rh tz - h t u t - r ji) e T] Z J dzdt = 0, 1 = 1,2, 
J (v 2 u ez C 2 + (u tt + p(h 2 ) z - f) e C) dzdt = 0. 



Now, we are allowed to take rj = phu and £ = w et in the above equalities. Hence 
taking in account that 

2 / [ hut hie dzdt = [ hf e (h,z)dz- [ h 2 e (0,z)dz, 1 = 1,2, 
Jo Jt Jt Jt 

2 / / u ett u e tdzdt= / u 2 t (t 1 ,z)dz - / u 2 t (Q,z)dz, 
Jo Jt Jt Jt 

2 / v 2 u tz u tzt dzdt = I v 2 \z) u 2 z (ti, z) dz — / v 2 {z) u 2 z (0, z) dz, 
Jo Jt Jt Jt 

and passing to the limit as e —> + , we obtain in analogy with (|4,27[) 

\ [ {ph 2 (t 1 ,z)+u 2 {t 1 ,z) + v 2 u 2 z {t 1 ,z))dz-^ f {phl + u 2 + v 2 u 2 z )dz 

* JT ^ JT 

= pr(j ■ h z — hi) dzdt + / / fu t dzdt 

Jo Jt Jo Jt 

<\ f 1 f pr{j 2 -h 2 z )dzdt+ (I \fu t \dzdt. 

z JO JT JO JT 

3. Finally, from P~4"T1) it follows that 



(4.47) 



u t (ti , z) dz < / (ph +u 1 + i> u 0z ) dz 
Jt 

+ [ [ prj 2 dzdt + 2 [ [ \fut\dzdt. 
Jo Jt Jo Jt 

Then, integrating the above inequality with respect to the time variable from 
to ti, we have 

nu 2 dzdt <2ti ( (phl + u\ + v 2 u 2 0z ) dz 
Jt 

+ 2tx [ [ prfdzdt + Atj f f f 2 dzdt, 
Jo Jt Jo Jt 

where we have used Cauchy's inequality. Again by Cauchy's inequality with 
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e = ti in (|4.47l) and the latter result, we obtain 

(ph 2 (t!) + u 2 t {ti) + y 2 u 2 z (h)) dz< (phi + u\ + v 2 u\ z ) dz 



.2. ,2,, ... - / ^1,2 „2 , „•>„- 

IT 

P r {f - dzdt + 2t x ( ( f 2 dzdt. 
io Jt 2 J J t 



□ 

Similarly we have done in the proof of the Existence theorem, item (7) , we 
obtain from (|4.45|) 

max \h\ + \\h\\ 2 ,u T + \\u\\ w i, HUt) <C 5 , heC a " a ^ 2 (U T ), (4.48) 

where C 5 = C 5 (T) > and a 2 € (0, a]. 

Now, let us show the uniqueness result. 

Theorem 4.6. (Uniqueness theorem). Under conditions of Theorem \4-l\ 
the periodically Cauchy diffraction problem (|3.1ip - (|3.16l) has a unique weak so- 
lution. 

Proof. 1. First, let h^,Uk, (k = 1,2), be two weak solutions of problem (|3.1ip - 
(|3.16p . well understood, for the same given initial data. For convenience, we 
define 

v := h,2 — hi and w := u 2 — U\. (4.49) 

Therefore, the functions v and w form a weak solution of the homogeneous 
problem 

v t = (rv z - h 2 w t - u u v) z in II T , (4.50) 

w t t = {v 2 w z -p(h 2 + hi) ■ v) z in II T , (4.51) 

v = 0, w = w t = on {0} x T. (4.52) 

Moreover, v and w satisfy the compatibility conditions (|3.15l) and (|3.16[) at the 
jump points. Then, we have the following equations satisfied 

( - vi r] t + (rvi z - hi2 w t - uu «j) rjz) dzdt = 0, 1 = 1,2, (4.53) 
( - w t Ct + v 2 w z ( z + p((h 2 + hi) ■ v) z () dzdt = 0. (4.54) 

2. We proceed as in the proof of Lemma T4. 5 1 Hence by analogy with (|4.46[) , 
we obtain from (|4"331 . (|4~5lj) 

T-e _ 

(vietV + (r v u - h i2 w t - uitvi) e f) z ) dzdt = 0, I = 1, 2, 

(4.55) 



JT 
T— 



/ / (v 2 w f _ z C, z + {w u +p{{h 2 + h 1 )-v) t C)dzdt = Q. 
Jo Jt 



22 



Moreover, for t\ 6 (0,T — e], we choose conveniently f), C, as 

m- 

and 



V ie (t), ifiG(0,ti] 
0, otherwise, 



C(t) 



w et (t), iftG(0,ti] 
0, otherwise, 



thus we observe that 

/ vutVie dzdt = / vf e (ti)dz- I vf e (0)dz, 1 = 1,2, 
o Jt Jt Jt 

2 / Wett Wet dzdt = / w%(ti)dz- / w^ t (0)dz, 
Jo jt Jt Jt 

/ w e tw e dzdt = I w^(ti)dz— / w^(0)dz. 
o Jt Jt Jt 

Now, we multiply the former equality in (|4.55l) by p, sum with the second one, 
and with such choice and observations, it follows that 

|ll^l)lll,T + ^lkt(tl)lli,T + ^lk^(*l)lll,T+ [J prv 2 z dzdt 
tl . 

P (uu vi z vi - u u v 2z v 2 - hnz V\ w t - h 21z v 2 w t ) dzdt, (4.56) 

o Jt 

where we have passed to the limit as e — > + and used the homogeneous initial 
data (|4.52l) . The result (|4.48p , obtained from Lemma l4~5l is now used to estimate 
the integral in the right hand side of (|4.56[) . Then, we have 

t_ _ 

V ("it V\z Vi + uu v 2z v 2 + hu vi w t + h 2z v 2 w t ) dzdt 

JT 

< pC 5 max|u|( ||wi z || 2 ,n T + ||v 2 z|| 2 ,n T + ||wt|| 2 ,n T 

LIT V 

Therefore, for almost all t\ € [0, T], we have from (|4.56l) and the above inequality 
2lK*0ll2,T+ 2lK(*i)ll2,T+ 211^(^)112^+ / / prv 2 z dzdt 



o Jt 



< 



C 5 max|u|( ||u z || 2 ,n T + IK|| 2 ,n T )- (4.57) 

Fly \ / 



3. Finally, let us show that v = and w = in Ht- Consider the equation 
(|4.50p . rewritten as 

v t + (u lt v) -(rv z ) z = (-h 2 w t ) 
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and recall the homogeneous initial condition v(0) — 0. Thus, if ||/i2^t||2,n T = 0' 
then the unique solution is v = 0. Consequently, by equation (|4.51| with 
w(0) = wt(0) = 0, we have w = 0, and we are done. On the other hand, if 
||^2W t ||2,n T = nj > 0, then from Theorem 8.1 [T31 P-192], again interpreting 
Q = [0,T] x [-1,2], Q' = [0,T] x [0, 1] and by the periodicity procedure, there 
exists a constant Cq(T ) > 0, which depends only on T whenever T £ (0,T ], 
such that 



maxK^I<a>(T ), VTe(0,T ]. 



Then, it follows that 

max|v| < ||fe2||2,n T ||wt||2,n r C 6 (T ) 
<C 5 (T )C 6 (T )\\w t \\ 2 ,n T , 
for each T € (0, To]. Hence from (|4.57[) and the above inequality, we obtain 



1 , 1 



ti 



p«(*l)lllT + dM*l)lli T + 5lK*(*l)ll2,T+ / prvidzdt 



< 



JT 



pC 5 2 (r )C 6 (T )(||« 2 || 2 ,n T ||^||2,n T + ||w t |l|n T ) 



Now, applying the Cauchy inequality with e = 2r (£5(^0)) 2 (C6(T )) , we 
obtain for almost all t x e [0, T],Te (0, T ] 



|ll«(ti)ll2,T + 5lK(*i)ll2,T + ~IK*(*i)ll2,T <PC 7 (T ) \\vH\\l,n T > 

where 



C 7 (T ) = Ci(T ) C 6 (T ) (l + C|(r0 4 ^ 6(r0) ) 
and integrating with respect to ii from to T, with T € (0,To], we have 

flMla,n T + ^IKH5,n T + y KII§,n T <^o C 7 (T ) IMlW ( 4 - 58 ) 
The inequality (|4.58[) is false when 

pT C 7 (T ) < 1/2, 
whenever wt, w z and v are non-zero functions. Consequently, taking 

Ti =min{(4pC7 7 (r )) _1 ,r }, 

we must have v = 0, w = in IL^ . The result follows applying the above 
procedure recursively, that is, after a finite number of steps, we get v = 0, w = 
in IT T . □ 
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4.3 Stability of weak solution 

The aim of this section is to study the stability of the weak solutions of the 
periodically Cauchy diffraction problem (|3.1ip - (|3.16p . with respect to variations 
of all coefficients, the free terms and the initial-data. 

In Sections 14.11 and 14.21 we have shown that the Cauchy problem (|3.11[) - 
(|3.16[) has a unique weak solution 

(heV 2 (ii T ), «6^ u (n r )), 

satisfying some additional properties, for instance 

max|Ai| + \\h z \\ 2 ,n T + IMIw^rnr) < C a , (4.59) 

where the positive constant Cs does not depend on h,u. Therefore, for each 
to e N we consider the following sequence of problems associated to (|3.11[) 

hf- = (r m h™ - h m ul n - r m j m ) in II T , (4.60) 

u m tt = {[y m f u™ - p{h m f) + f m in n T , (4.61) 

h m = h™ on {0} x T, (4.62) 

u m = uT = «y on {0} x T. (4.63) 

In fact, by a standard density argument, it is enough to suppose that 

r" l (z), j m (t,z), h%(z), vZ(z), u?(z), f m (t,z) 

are smooth functions satisfying the conditions of the uniqueness and existence 
theorems. Then, for each mgff, there exists a unique weak solution 

(h m eV 2 (lL T ), u m e ^ 2 M (n T )), 

satisfying 

max|^| + \\h™h, UT + IKIUm (IIt ) < C 9 , (4.64) 

LIT 2 

where the positive constant Cg does not depend on h m , u m and m. More- 
over, the transmission conditions (I3.15[) — (|3.16|) can be dropped owing to the 
smoothness of the solution. 

Theorem 4.7. (Stability theorem). Suppose that the sequences {r m } and 
{v m } are uniformly bounded and converge a.e. to r,v respectively. Also the 
sequences {j m }, {/ m }, {h™}, {u™}, {u™} converge to j , / , Hq,uq,ui in the re- 
spectively norms of the spaces to which they belong according to conditions of 
Theorem \4-l\ Then, the sequence of weak solutions {(h m , u m )}^ =1 , 

(h m G V 2 lfl {U T ),u m G VK 2 14 (n T )) (Vm > 1), 
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of the associated problems (|4.60p - (|4.63p converges in such spaces to the weak 
solution (h,u) of the limit problem (|3 . 11 [) — f|3 . 1 6|) . 

Proof. We argue by contradiction. Therefore, let us assume that the sequence 
{(h m , M m )}m=i does not converge to (h,u), and since it is uniformly bounded, 
we can extract a subsequence, namely {{h mi ,u m " i ^}°^ 1 , such that 

h m * h weakly in £^(0, T; L 2 (T)), (4.65) 

h™* h z weakly in L 2 (U T ), (4.66) 

u m * -v {j weakly in 1^(0, T; W^(T)), (4.67) 

u m * w< weakly in £^(0, T; L 2 (T)). (4.68) 

By the estimates obtained from the existence and uniqueness theorems, it is not 
difficult to show that, for each 77, £ £ W 2 ' (Ht) 

YC"" 1 rjt dzdt — > — / / hife (izcft, 



JJut 



r h71 i r\ z dzdt — > // rhi z rj z dzdt, 1=1,2, 



n 3 



(i/™*) 2 1*^** C« -> // v 2 u z ( z dzdt. 

Moreover, if we proceed as in the proof of the existence theorem, then we can 
show that 

/i™ 1 it™* ?7 Z dz<ii > — / / hiurj z dzdt, 1 — 1,2, 
n T JJn T 

p((h m ') 2 ) z Cdzdt^ ff p(h) z (dzdt, 

for all rj, ( € W2' 1 (Ht) ■ Consequently, from the convergence above we conclude 
that (h, u) is a weak solution weak the Cauchy problem (|3.11l) - (|3.16[) with initial 
data ho, uq and u\. But, by uniqueness of the solution, we have 

h = h and u = u, 

which is a contradiction with our initial assumption. □ 



Besides we have proven the stability result, we give a more refined result con- 
cerning the estimative. The result is established in the case when v is supposed 
to be a smooth enough function. 
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Theorem 4.8. Let h,u and h m ,u m , for each m S N, be weak solutions of 
problems (| 3 . 1 1 1) (| 3 . 1 6 [) and the associated (|4. 60[) — (|4.63[) respectively, and define 

v m := h m - h, w m := u m - u. 

Then, there exists a positive number 5, independent of h,u,h m ,u m , ij.,i2, such 
that for any t\,t2 € [0,T], < t<i — t\ < S, the following inequality holds true 

vrai max {p |K||| T + ||<*||2, T + 2 ||^<||| T } + 2pr ||<|||n tlt2 
te[ti,t 2 ] I- - 1 12 

< P C 10 {||(r m - r)^||| ntit2 + \\r m j m - r ,n tlt2 } + 2p ||«"*(-, *i)||i, T 

+ II/'" - f\\l,n Ht2 + 2 |K(*i)||!,t + 2 \\vw?(h)\\l r , (4.69) 
where n tlt2 = (iijig) x T, and Cio is a positive constant independent oft±,t2- 
Theorem 14.81 is proved in much the same way as Lemma 4.2 from |20] . 
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